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ABSTRACT 


This  project  developed  and  evaluated  a  new  mathematical  decompression  model  with 
asymmetrical  gas  kinetics.  The  intended  application  was  to  support  U.S.  Navy  diving  operations 
involving  surface  decompression  with  oxygen  (O2)  following  air  diving  (Sur-D  O2).  Before  this 
effort,  the  most  sophisticated  model  for  predicting  human  risk  of  decompression  sickness  (DCS) 
following  exposure  to  elevated  O2  mixtures  was  one  developed  in  1998  and  based  on  only  a  small 
amount  of  Sur-D  02  data.  The  present  effort  added  more  than  4,000  dives,  particularly  dives  with 
high  02  and/or  Sur-D  O2,  to  the  data  for  calibration  of  the  new  model.  About  half  of  the  added  dives 
were  experimental  exposures  involving  rats,  in  the  hope  that  higher-risk  animal  dives  would 
improve  prediction  accuracy  for  higher-risk  human  dives.  It  was  also  thought  that  sharing 
parameters  between  the  species,  particularly  those  parameters  defining  the  effect  of  O2  on 
decompression  risk,  might  enhance  model  performance.  However,  we  were  unable  to  demonstrate 
an  advantage  of  the  rat-human  model  over  the  human-only  or  the  1998  models.  We  used  all  three 
models  to  evaluate  three  possible  alternative  Sur-D  02  procedures,  each  alternative  having  more 
flexibility  than  those  currently  in  use  but  varying  in  amounts  of  DCS  risk.  The  Navy  can  now 
evaluate  changes  in  Sur-D  62  procedures  by  using  all  three  models.  When  new  procedures  are 
considered  for  Fleet  use,  this  approach  may  enhance  decision  making  without  requiring  a  manned 
dive  trial. 
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INTRODUCTION 


Recent  experience  shows  that  the  Navy  still  lacks  the  ability  to  effectively  employ  02  as 
a  decompression  tool.  This  is  evident  by: 

1 .  perceived  high  decompression  sickness  (DCS)  incidence  using  surface 
decompression  procedures  with  02  following  air  diving  (Sur-D  02)  during  the 
recovery  of  TWA  flight  800, 1  and 

2.  NAVSEA’s  acknowledged  continued  desire  to  improve  the  use  of  02  in  order  to 
optimize  decompression  schedules  ( balance  in-water  vs  total  decompression  time 
depending  upon  objective,  etc). 

In  order  to  evaluate  changes  to  decompression  tables,  tools  are  needed  that  can  predict 
DCS  outcomes  based  upon  the  decompression  schedule.  For  the  past  two  decades, 
probabilistic  models  have  been  developed  for  this  purpose.2,3  Until  now,  the  most 
sophisticated  probabilistic  model  used  to  evaluate  procedures  using  elevated  partial 
pressures  of  02  (including  Sur-D  02)  was  based  on  only  a  small  amount  of  Sur-D  02 
data.3  A  new  probabilistic  model  was  developed  in  present  work  and  used  to  predict 
human  DCS  risk  following  exposure  to  elevated  02  mixtures.  The  three  notable  aspects 
of  this  model  are 

1 .  time  constant  asymmetry  between  uptake  and  washout, 

2.  an  expanded  human  nitrox  data  set,  and 

3.  the  use  of  animal  data  to  expand  depth/time/P02  scope  of  calibration  data 
beyond  that  of  the  human  data. 


METHODS 


MODELS 

We  used  a  class  of  models  that  incorporate  asymmetrical  gas  kinetics  and  that  this 
report  refers  to  as  “ASYM.”  Based  on  standard  exponential  kinetics  for  the  gases,  the 
ASYM  model  was  the  basis  of  earlier  work  with  human4  and  human/rat  decompression 
models.5  The  asymmetry  has  been  used  previously  in  a  deterministic  model4,  but  this  is 
a  first  attempt  to  estimate  the  parameters  in  a  probabilistic  context.  We  chose  this 
model  over  the  Linear-Exponential  (LE)  models  of  previous  work  due  to  its  flexibility  in 
allowing  gas  elimination  to  be  either  faster  or  slower  than  uptake,  which  was  not 
possible  for  the  LE  models. 

Our  models  were  constructed  with  a  series  of  theoretical  gas-exchange  “compartments” 
or  “tissues”  which,  for  modeling  purposes,  were  independent,  parallel-perfused 
compartments.  The  number  of  these  compartments  affects  how  well  the  model  fits  the 
data.  The  novel  feature  of  ASYM  is  that  the  time  constant  governing  blood-tissue  gas 
exchange  kinetics  for  compartment  j  assumes  one  of  two  possible  values  depending  on 
the  calculated  tissue  inert  gas  overpressure  in  that  compartment: 
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(1b) 


(la) 


Here  tau  (t)  is  the  compartment  (/)  time  constant  of  gas  7,  P,  is  the  partial  pressures  of 
the  inspired  gas  7’,  Pamb  is  the  ambient  pressure,  P thr  is  the  estimated  threshold 
parameter,  and  asym  is  the  asymmetry  parameter  affecting  the  rate  of  gas  washout. 
The  modified  time  constants  due  to  the  asym  parameter  are  active  during 
decompression;  thus,  this  parameter  can  be  thought  of  as  either  slowing  down  ( asym  > 
1 )  or  speeding  up  ( asym  <  1 )  gas  washout  during  this  phase  of  a  dive.  With  asym  =  1 , 
gas  uptake  and  washout  are  equal. 

This  project  was  aimed  at  decompressions  using  high  levels  of  02,  so  an  explicit 
contribution  by  02  to  DCS  risk  was  incorporated  into  the  model.  This  was  done  in  a 
manner  similar  to  that  of  Parker  et  al.3  for  human  diving  in  their  model  #2.  That  model, 
which  will  be  referred  to  as  “JAP98,”  defined  the  role  of  02as  adding  to  the  inert  gas 
load  by  contributing  to  bubble  formation  and/or  growth  when  the  inspired  02  partial 
pressure  is  greater  than  parameter  P set-  Thus,  the  portion  of  the  ambient  partial 
pressure  of  02  that  exceeds  the  Pse(  value  is  treated  as  an  inert  gas.  In  our  model,  the 
gas  exchange  between  the  ambient  environment  and  the  tissue  was  governed  by  blood 
perfusion  limitations  according  to  the  equation 


(2) 


where  P(.  is  the  partial  pressure  of  gas  /  dissolved  in  the  tissue,  ?amb  is  the  partial 

pressure  of  gas  /  in  the  inspired  gas,  and  r,  is  the  time  constant  for  gas  i,  all  in 
compartment  j.  For  the  case  of  02,  the  Yamb  is  replaced  by  (P  .  -  P set),  with  this  term 

constrained  to  be  greater  than  zero: 


dP, 


(3) 


To  evaluate  the  ASYM  model,  we  needed  to  define  risk.  The  instantaneous  risk  for  a 
given  compartment  was  defined  to  be  the  supersaturation  above  a  threshold  (P thr): 


(4) 


where  gain  is  a  scale  factor  and  j  is  the  index  of  the  individual  compartments.  The 
probability  of  no-DCS  in  the  0  -  T  interval  is  then 
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(5) 


P(noDCS)=e 


1 

|  risk  »dr 
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where  integration  of  the  risk  occurs  until  risk  falls  to  zero.  We  included  the  onset  time  of 
DCS  in  our  models,  as  experience  suggests  that  including  this  time  improves  estimates 
of  model  parameters.6  Similarly,  the  probability  of  DCS  occurring  within  a  time  interval 
between  times  Ti  and  T2  is  defined  as 


P(DCS)-- 


Tl 

- 1  risk  •  dt 


f  ^  \ 

-  risk •  dt 

l-e  Tl 


V 


/ 


(6) 


where  Ti  is  the  last  time  the  individual  was  known  to  be  symptom  free  and  T2  is  the  first 
time  that  symptoms  are  sufficient  to  make  a  diagnosis.6 

Once  we  have  defined  the  probability  of  an  event  occurring,  we  can  use  the  method  of 
maximum  likelihood  to  match  this  to  the  outcomes  for  these  profiles,  which  assumes 
independent  events.  The  likelihood  (L)  is  the  joint  probability  of  observed  outcomes 
given  the  model7: 


L  =  flP(DCS„y"  •P(noDCSjv'\,  (7) 

n 

which  is  the  product  of  probabilities  for  each  dive  profile  ,  n,  over  all  the  dive  profiles,  N, 
where  y/n  =  1  if  DCS,  or  if/n  =  0  if  no  DCS.  As  L  becomes  a  small  number  when  multiple 

dives  are  considered,  we  used  the  log-likelihood  (LL),  which  is  the  natural  logarithm  of 
the  above  expression.  The  best  match  between  a  model  and  the  outcomes  in  the  data 
occurs  when  the  likelihood  is  maximized.  We  used  parameter  estimation  techniques  to 
modify  the  model  parameters  to  find  this  maximum  value,  which  is  our  “best  fit.” 


DATA 

Human 


A  total  of  6,207  human  nitrox  exposures  were  used  for  (Table  1 ),  nearly  2,000  more 
than  in  any  prior  published  analysis.  These  additional  exposures  included  a  number 
involving  high  inspired  partial  pressures  of  O2  that  have  never  been  included  in  a 
calibration  data  set  for  a  probabilistic  model.  Specifically,  the  data  set  consisted  of 

1 .  the  26  data  sets  used  with  a  previously  reported  human  model  that  included  O2 
in  the  prediction  of  DCS,3  (this  information  is  shown  on  Table  1 ,  [3]),  and 

2.  an  additional  12  data  subsets  from  earlier  research  efforts.8 
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Table  1.  Human  data  sets  used  for  modeling[1]. 


Data  Set 

Depth 

(fsw) 

Bottom  Time 
(min) 

Number  of 
Exposures 

Observed  Cases 
of  DCSP1 

Single  Air 

EDU885A131 

50-190 

14-244 

483 

30 

UPS290 

60  -  625 

0-5 

299 

4 

DC4Wl3J 

49  -  265 

9-100 

244 

8.4 

EDU849LT2l4J 

100-150 

27-60 

73 

9.2 

NMR97NOD 

100-170 

38-85 

103 

3.4 

NMRNSWl3J 

61.5 

82  -  364 

91 

5.5 

PASA131 

101  -  151 

30-62 

72 

5.2 

SUBX87l3J 

80  -  602 

0-2 

58 

2 

Single  Non-Air 

NMR8697T3r 

25-130 

30  -  240 

477 

12.8 

EDU1 180Sra 

75-150 

38-126 

120 

10 

EDU885M13J 

100-150 

44-66 

81 

4 

Repetitive  and  Multilevel  Air 

pamla[3J 

51  -82 

61  -452 

236 

14.2 

EDU885ARL3J 

80-150 

29-62 

182 

11 

PARA131 

61  -  152 

25-70 

135 

7.3 

DC4WRl3J 

120-180 

30-40 

12 

3 

Repetitive  and  Multilevel  Non-Air 

EDU184131 

40-150 

23  -  366 

239 

11 

PAMLOAS131 

61  -82 

33  -  637 

140 

5.3  1 

PAMLAOD131 

61  -82 

241  -  540 

134 

6 

EDU885S131 

60-150 

43-124 

94 

4 

Air  with  Oxygen  Decompression 

NMR94EODl3J 

60-170 

40  -  240 

284 

17.9 

DC8AOD13J 

59  -  300 

2-60 

256 

3.2 

DC8AOW131 

90-180 

2-208 

46 

3.1 

[1)  Listed  by  original  data  set  names 

(2]  Observed  cases  =  diagnosed  cases  +  0.1  •  marginal  cases 

(31  Data  sets  to  which  model  JAP98  was  fit,  as  described  in  Parker  et  al.,  1998  (reference  3). 
[4'  Only  the  resting  dives  were  used  from  this  study. 
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Table  1.  Human  data  sets  used  for  modeling  (continued). 


Data  Set 

Depth 

(fsw) 

Bottom  Time 
(min) 

Number  of 
Exposures 

Observed 
Cases  of  DCS121 

Saturation 

ASATNSMlJJ 

25-132 

2880-  11760 

132 

20.1 

ASATEDU131 

50-60 

740  -  2884 

120 

15.7 

SUREX 

65-75 

5460  -  7020 

24 

5.3 

ASATNMRiaj 

20-24 

4320-6181 

50 

1 

NMR9209 

20-22 

4284  -  4399 

48 

2.5 

EDUAS45 

33-99 

360-2160 

26 

5.6 

ASATDC 

26-33 

1440-1454 

23 

8.1 

ASATARElJ| 

23-78 

2862  -  2895 

165 

21.3 

Subsaturation 

EDU849S2l4J 

33-40 

720 

30 

5.8 

NSMeHR131- 

28-40 

359  -  360 

57 

3.2 

Surface  Decompression  on  Air 

EDU545SUR141 

100-170 

30-85 

102 

10.3 

NMRASUR90 

60-120 

49-180 

64 

0 

Surface  Decompression  on  C 

•xygen 

EDU1351SDlbJ 

74-214 

20-130 

1035 

43.3 

DC8ASUR181 

59  -  236 

30-70 

358 

10.1 

DCSUREP181 

147-172 

30  -  230 

69 

1 

NMROSUR90 

60-120 

49-179 

45 

1 

[3]  Data  sets  to  which  model  JAP98  was  fit,  as  described  in  Parker  et  al.,  1998  (reference  3). 

[4)  Only  the  resting  dives  were  used  from  this  study. 

151  Sur-D  02  was  assumed  to  have  inspired  02  fraction  of  90%. 

For  previous  efforts,  we  had  established  a  strict  set  of  requirements9  for  acceptance  of 
calibration  data.  One  of  these  requirements  was  that  we  reconcile  data  questions  with 
the  author/principal  investigator  of  the  study  that  generated  the  original  data  set. 
However,  the  historical  data  was  too  old  for  such  reconciliation.  Therefore,  we  decided 
to  relax  this  specific  criterion.  We  were  confident  that  the  depth-time-gas  mix  histories 
could  be  reconstructed,  but  were  less  certain  about  addressing  questions  on  DCS 
outcomes.  The  severity  of  the  historic  dive  profiles  resulted  in  high  DCS  incidences17 
that  would  probably  be  even  higher  than  recorded  if  today's  more  conservative 
diagnostic  criteria  had  been  used.  Fortunately  the  postdive  descriptions  of  DCS  signs 
and  symptoms  were  available,  so  a  panel  of  DMOs  could  (and  did)  rediagnose  all 
outcomes  with  a  standard  set  of  criteria.8,9 

The  human  data  showed  an  average  DCS  incidence  of  4.8%  (1 1 .2%  when  marginal 
cases  were  included  as  full  cases),  with  a  depth  range  of  20  to  625  fsw  (maximum  of 
297  fsw  when  submarine  escape  dives  are  excluded),  and  an  O2  partial  pressure  range 
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of  0.3  to  4.2  ata  (max  2.7  ata  with  the  exclusion  of  the  submarine  escape  dives,  which 
are  characterized  by  very  short  duration  to  great  depths,  with  rapid  ascents). 


Rat 


A  total  of  2,390  rat  dives  used  for  the  animal  component  included  many  dives  with 
elevated  O2: 

1 .  2,120  nitrox  dives  over  a  large  range  in  pressures,  bottom  times,  decompression 
profiles,  and  02  pressures,  and 

2.  270  dives  from  an  experiment  that  involved  switching  the  inspired  partial 
pressure  of  02  at  depth  before  decompression. 

The  2,120  nitrox  dives  had  a  DCS  incidence  of  50.5%,  with  a  depth  range  of  141-275 
fsw  and  oxygen  partial  pressures  up  to  3  ata.  The  270  oxygen  switching  dives  had  a 
61 .1%  incidence  of  DCS,  with  a  depth  range  of  225  to  275  fsw.  These  dives  were 
conducted  for  60  min  on  air,  with  135  control  animals  continuing  on  air,  while  the  others 
switched  to  42%  nitrox  before  decompression. 

The  2,120  nitrox  dives  have  been  described  previously.10'14  The  270  O2  switching  dives 
were  conducted  in  a  similar  fashion. 

Time  of  DCS 


The  12  new  human  data  sets  were  found  to  contain  inconsistencies  in  the  times  related 
to  onset  of  DCS.  To  correct  this  problem,  we  chose  a  conservative  approach  for 
standardizing  the  Ti  and  T2  values  for  the  human  dives.  The  Ti  values  were  set  to  an 
earlier  time  than  recorded:  either  the  closest  documented  medical  check  time  or  the  end 
of  the  bottom  time,  whichever  was  later.  The  T2  times  were  moved  to  a  later  time  than 
recorded,  which  was  the  next  documented  medical  check  time.  As  our  rat  data 
contained  only  the  time  at  which  DCS  was  first  observed  during  the  experiment,  we 
needed  to  establish  rules  for  defining  T 1  and  T2  times  for  these  dives.  As  with  the 
human  data,  we  were  conservative  with  the  rat  data:  we  defined  Ti  for  rats  as  the  time 
leaving  the  bottom  and  T2  as  the  original  time  when  DCS  had  been  observed. 

These  modifications  reduced  the  effect  of  the  additional  dives  on  the  estimated 
parameters  in  the  model,  as  they  made  the  times  less  precise.  However,  without 
evaluating  and  setting  the  T  times  according  to  a  fixed  and  consistent  set  of  rules  with 
the  full  rigor  that  was  used  for  the  previous  data  sets,  we  felt  that  it  was  necessary  to  err 
on  the  side  of  caution. 
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PARAMETER  ESTIMATION 


Model  parameter  values  were  adjusted  to  maximize  the  model  log  likelihood  (LL)  about 
the  calibration  data  using  a  modified  Marquardt  nonlinear  estimation  algorithm.  This 
fitting  process  required  defining  the  starting  values  of  all  parameters.  However,  for  such 
complex  models  as  these,  we  had  to  perform  the  fitting  process  repeatedly  with  many 
possible  combinations  of  starting  parameter  values  to  improve  the  chances  of  finding  a 
global  maximum  in  the  likelihood  surface  (i.e. ,  best  fit).  To  facilitate  the  fitting  process, 
we  created  a  tool  that  generated  multiple  sets  of  starting  parameter  values  randomly 
selected  from  a  user-defined  range.  This  tool  helped  us  perform  a  quasi-global  search 
of  the  parameter  space  with  the  starting  parameter  values  limited  to  plausible  areas. 
Once  we  obtained  a  possible  maximum  in  the  likelihood  surface,  we  had  to  test  whether 
it  was  truly  a  global  or  a  merely  a  local  maximum.  This  was  done  by  modifying  the  last 
best-fit  parameter  values  by  randomly  adding  or  subtracting  varying  amounts  to  create 
new  starting  parameter  value  sets,  which  were  then  used  to  conduct  confirmation 
searches  for  other  maxima.  If  none  were  found,  we  broadened  our  search  until  we  felt 
confident  that  the  best  fit  was  obtained.  More  than  500  runs  were  completed  to  confirm 
each  of  the  final  parameter  sets  reported. 

For  the  combined  species  model,  each  rat  parameter  was  set  equal  to  the 
corresponding  human  parameter  plus  a  delta  parameter.  For  example,  the  rat  asym 
parameter  was  given  by: 


asymra,  =  asym \human  +  A  asymral  (8) 

where  the  possible  estimated  parameters  were  asymhuman  and  A asymral .  In  order  to  link 

the  species  in  the  combined  species  model,  we  defined  several  of  the  parameters  as 
common  (i.e.,  set  equal)  between  humans  and  rats  by  fixing  the  delta  parameter  at 
zero.  This  minimized  the  number  of  estimated  parameters,  and  provided  a  means  by 
which  the  rat  data  could  affect  model  fit  to  the  human  data.  These  shared  parameters 
were 

1 .  asym ,  the  asymmetry  parameter  which  affects  the  mathematical  relationship 
between  gas  washout  and  gas  uptake,  and 

2.  Psef,  which  defines  the  pressure  at  which  O2  begins  to  contribute  to  risk. 

The  parameters  for  which  the  deltas  were  estimated  were: 

1 .  r,  time  constants  for  N2  and  02, 

2.  P thr,  the  threshold  pressure  above  which  risk  accumulates,  and 

3.  gain,  the  scale  factor  for  each  compartment. 

Estimation  of  these  deltas  effectively  made  these  parameters  species-specific. 
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The  likelihood  ratio  ( LR )  test  was  used  to  evaluate  the  significance  of  estimated 
parameters  on  the  basis  of  improvement  in  fit.7  As  parameters  were  added  to  the 
model,  they  were  tested  for  significant  improvement,  by  calculating  the  LR  statistic 

LR  =  2  •  Km, ££*»].  (9) 

where  LLbase  is  the  LL  of  the  model  prior  to  the  parameter  additions,  and  LLadditionai  and  is 
the  LL  of  the  same  model  with  the  additional  parameters  being  estimated.  This  setup 
creates  a  nesting  relationship  for  the  models,  making  one  a  subset  of  the  other.  This 
allows  us  to  look  up  the  LR  statistic  on  the  Chi-squared  tables  (with  the  number  of 
additional  parameters  being  the  degrees  of  freedom)  to  assess  if  the  addition  is 
statistically  significant  at  a  given  confidence  level  (we  use  95%).  In  our  case,  all  of  the 
parameters  exist  in  the  model  code  at  the  start,  but  had  their  values  fixed  to  remove  any 
contribution  to  the  model  risk  (thresholds  fixed  at  zero;  PsetS  fixed  at  99). 

As  a  final  step  in  the  addition  of  parameters,  we  tested  the  addition  of  delta  asym  and 
Pse(  parameters  for  rats,  which  effectively  separated  the  species  by  allowing  all  possible 
parameters  to  have  independent  values.  The  LR  test  of  the  addition  of  these  two 
parameters  is  effectively  also  a  test  that  determines  whether  the  human  and  rat  data 
sets  were  statistically  combinable: 


LR  =  2  • 


[LL 


independen  t 


LL 


combined 


] 


(10) 


where  LLcombjned  is  the  LL  of  the  described  combined  model  fitted  to  the  combined  data 
sets,  and  LLjndependent  and  is  the  LL  of  the  same  model  with  remaining  delta  parameters 
estimated,  which  effectively  made  the  human  and  rat  expressions  independent.  Data 
sets  not  combinable  by  this  criterion  generally  produce  better  predictions  of 
decompression  risk  when  modeled  separately.  However,  past  work,  with  different  types 
of  human  dive  trials,  has  succeeded  while  ignoring  this  fault.9 


RESULTS 


MODEL  PARAMETERS 

Model  parameters  estimated  for  the  human-only  and  combined  human/rat  data  are 
presented  in  Tables  2  and  4.  Only  parameters  found  to  be  significant  at  P<  0.05  by  the 
LR  test  are  reported.  Parameter  significance  was  determined  by  building  on  the  basic 
model  (with  minimal  parameters),  adding  one  or  two  parameters  at  a  time,  and  then 
using  the  appropriate  number  of  degrees  of  freedom  (one  per  parameter  added)  to 
apply  the  LR  test.  The  02  parameters  (Pse(  andr)  were  added  to  determine  how 
significant  any  02  effect  would  be,  as  demonstrated  by  a  significant  increase  in  LL. 
Because  of  the  sequential  nature  of  the  model  fitting  process,  the  true  confidence  level 
of  both  the  parameters  and  the  final  model  should  be  less  than  what  might  be  expected. 
In  some  cases,  parameters  were  fixed  at  specific  values  rather  than  estimated,  and 
these  parameters  are  noted  by  having  no  standard  error. 
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The  best-fit  human-only  model  (Table  2)  had  three  compartments  with  no  asymmetry 
and  a  significant  02  effect  only  in  compartment  1 ,  with  10  significant  estimated 
parameters  out  of  a  possible  18.  The  human-only  model  is  a  presentation  of  just  the 
human  parameters  of  the  independent  model,  as  it  is  equivalent  to  fitting  just  the  human 
parameters  to  the  human-only  data,  for  the  purposes  of  making  risk  predictions  on 
human  exposures.  Because  the  estimated  value  for  Pse(  was  small  in  compartment  1 
with  a  large  confidence  interval,  it  was  fixed  to  zero.  Fixing  the  Psej  value  improved  the 
ability  to  estimate  the  other  parameter  values  and  effectively  had  us  estimating  1 1 
parameters.  This  effectively  treated  02  in  a  fashion  similar  to  N2  in  the  definition  of  DCS 
risk  for  that  compartment.  To  remove  the  influence  of  02,  Pset  was  fixed  at  99.0  in  the 
other  two  compartments,  because  the  02  effect  was  not  found  to  be  significant.  This 
finding,  in  turn,  made  the  02  time  constants  meaningless,  and  values  in  these 
compartments  are  therefore  denoted  as  N/A.  Having  similar  values  for  the  N2  time 
constants  for  compartment  2  and  3,  with  rather  large  errors  that  allow  the  parameter 
ranges  to  overlap,  was  disconcerting.  However,  the  thresholds  for  these  two 
compartments  were  very  different  and  standard  errors  on  the  thresholds  were  tight. 

This  parameter  set  was  confirmed  by  more  than  1 ,000  additional  model  fits. 


Table  2.  Human-only  ASYM  model:  Parameter  values  +  standard  errors. 


Parameter 

Compartment 

Compartment 

Compartment 

1 

2 

3 

N2  Time  Constant  (min)  rN 

3.8±0.5 

330  +  20  . 

390 1  70 

02  Time  Constant  (min)  r02 

1.410.2 

N/A 

N/A 

asym 

1.0* 

1.0* 

1.0* 

Natural  Logarithm  of  the  Gain 

-5.410.4 

-7.510.1 

-4.710.3 

P thr  (atm) 

0.710.2 

0.07  ±  0.01 

0.51  1  0.04 

P set  (atm) 

0.0* 

99.0* 

_ i 

99.0* 

*  Fixed 


The  human  and  rat  data  sets  failed  the  test  for  combinability,  as  evidenced  by  the  very 
significant  increase  in  LL  (+  17.1)  with  the  addition  of  5  parameters  (2  APsets,  3  Aasyms) 
when  the  model  made  the  predictions  for  the  2  species  independent  of  the  calibration 
data  for  the  other  species  (Table  3).  However,  our  experience  has  shown  that  it  is  often 
difficult  to  pass  the  LR  test  when  assessing  the  combination  of  human  data  sets.  Given 
that  one  of  the  primary  goals  of  this  project  was  to  produce  and  evaluate  a  multispecies 
model,  we  proceeded  with  the  combined  human/rat  model. 
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Table  3.  Log  Likelihood  values  for  the  combined 
and  independent  species  ASYM  models. 


Data 

Log 

Likelihood 

(LL) 

Number  of 
Parameters 

Combined 

-2875.7* 

27 

Independent 

-2858.6 

32 

Human 

-1538.4 

11 

*,  Significantly  different  (P  <  0.01 )  from  separate  species  models 

The  best-fit  human/rat  model  (Table  4)  was  a  3-compartment  model  with  asymmetry  in 
all  compartments  and  an  O2  effect  in  compartments  1  and  2.  Using  the  likelihood  ratio 
test,  we  found  27  of  the  36  possible  parameters  in  the  combined  model  to  be 
statistically  significant.  The  observed  asymmetry  with  faster  washout  than  uptake  is  a 
phenomenon  observed  previously  in  rats.14  For  this  combined  model,  the  asym  and  Pse( 
parameters  were  forced  to  be  the  same  for  both  the  rat  and  the  human  compartments. 

Table  4.  Combined  human/rat  ASYM  model:  Parameter  values  +  standard  errors. 


Parameter 

Compartment  1 

Compartment  2 

Compartment  3 

Human 

Rat 

Human 

Rat 

Human 

Rat 

N2  Time 
Constant 
(min)  rA,2 

4.8  ±1.1 

0.6  ±0.1 

360  ±  30 

0.4  ±0.1 

620  ±  80 

10.6  ±0.3 

O2  Time 
Constant 
(min)  r0t 

2.4  ±  1.5 

80±  10 

0.9  ±0.4 

4  ±  2 

N/A 

N/A 

asym 

0.39  ±0.07 

Human** 

0.82  ±  0.07 

Human** 

0.32  ±0.04 

Human** 

Natural 
Logarithm  of 
the  gain 

-4.7  ±0.4 

0.5  ±0.3 

-7.4  ±0.1 

2.7  ±0.3 

-4.4  ±0.3 

-0.5  ±0.2  . 

Pf hr  (atm) 

0.4  ±0.2 

1.8  ±0.2 

0.05  ±0.01 

0.4  ±  1.1 

0.44  ±0.03 

2.9  ±0.1  • 

P set  (atm) 

0.2  ±0.1 

Human** 

1.1  ±0.2 

Human** 

99.0* 

Human** 

*  Fixed 

**  Fixed  equal  to  the  human  value. 
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The  values  of  the  human  parameters  for  the  models  fitted  to  the  human-only  and 
combined  data  were  different,  but  the  predictions  for  the  human  data  sets  reported  in 
reference  8  were  not  greatly  affected  (see  Appendix  A).  Thus,  including  the  higher-risk 
rat  dives  in  the  calibration  data  did  not  improve  the  ability  to  predict  human  dives.  That 
the  human-only  model  did  not  include  asymmetry  but  predicted  comparably  to  the 
combined  species  model  indicates  that  the  human  data  was  described  equally  well  with 
or  without  asymmetry. 

PREDICTED  DCS  RISKS  OF  PRESENT  AIR  AND  SUR-D  02  TABLES 

The  predicted  risk  of  DCS  was  not  observed  to  be  constant  over  the  entire  set  of 
present  Air  and  Sur-D  02  tables,  with  risk  as  low  as  2%  for  short  dives  and  >  5%  for 
longer  dives  (see  Appendix  B).  We  note,  however,  that  the  model  predictions  for  any 
schedule  are  for  the  worst  case  —  for  the  maximum  time/depth  —  as  we  cannot  foresee 
the  actual  depths  and  times  that  will  be  used.  Differences  in  risk  predictions  between 
in-water  and  surface  decompression  schedules  were  inconsistent,  although  some  of  the 
inconsistencies  may  have  resulted  from  the  different  techniques  used  to  develop  the 
two  sets  of  tables  (reference  16  (in-water)  and  reference  17  (Sur-D  02)). 

DISCUSSION 


With  an  emphasis  on  dives  with  high  02  and/or  Sur-D  02,  this  effort  added  over  4,000 
additional  dives  for  estimating  the  new  model.  About  half  of  these  added  dives  were 
experimental  exposures  involving  rats,  in  the  hope  that  the  higher-risk  animal  dives 
would  improve  prediction  accuracy  for  higher-risk  human  dives.  We  also  thought  that 
sharing  of  parameters  between  the  species,  particularly  those  parameters  defining  the 
effect  of  oxygen  on  decompression  risk,  might  enhance  model  performance. 
Unfortunately,  we  were  unable  to  demonstrate  any  significant  advantage  of  the 
multispecies  model  over  the  human-only  or  the  previously  reported  JAP98  model.3 

It  is  interesting  to  note  that  in  the  combined  human/rat  model,  the  asymmetry 
parameters  are  <  1  for  all  three  compartments  indicating  that  modeled  gas  exchange  is 
faster  for  washout  than  for  washin.  This  result  is  the  opposite  from  that  previously 
reported  for  linear/exponential-based  models  using  a  subset  of  the  present  human  data 
(see  Table  1 ,  JAP98).  Those  models  displayed  slower  off-gassing  characteristics  by 
using  bubble-like  dynamics.4  However,  there  is  some  experimental18  and  theoretical19 
basis  for  faster  elimination.  This  effect  is  driven  by  the  rat  data,  as  the  human  only 
model  did  not  support  assymetrical  gas  kinetics  (See  Table  3). 

We  did  not  recalibrate  the  JAP98  model  with  the  expanded  human  data  set.  However, 
both  the  JAP98  and  ASYM  models  agreed  well  when  comparing  the  model  predictions 
and  observed  outcomes  for  previous  data  sets  as  seen  in  Table  5.  Some  data  sets 
were  better  described  by  JAP98;  others  by  ASYM.  No  consistent  pattern  was 
discernible.  It  should  be  noted,  that  while  the  older  data  (EDU557,  EDU1351SD)  was 
rediagnosed  with  a  standard  set  of  criteria,  it  is  highly  likely  that  symptoms  of  DCS  that 
would  be  treated  today  (under  present  criteria),  would  not  have  been  noted  at  the  time; 
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thus  the  observed  incidences  for  these  data  sets  are  expected  to  be  low  compared  to 
the  models  predictions.  A  complete  comparison  of  the  data  sets  documented  in 
reference  8  appears  in  Appendix  A. 


Table  5.  Predicted  DCS  risks  for  previous  data  sets. 


Dive  Trial 

#  Dives 

#  DCS 

#DCS  Predicted  by  Model 
(95%  Confidence  Interval) 

JAP98 

ASYM 

Human/Rat 

ASYM 

Human-only 

1957  Standard  Air 

568 

27 

36 

26 

28 

(EDU557) 

(29  -43) 

(22  -  30) 

(23-32) 

1951  Sur-D  02 

1035 

43 

66 

52 

55 

(EDU1351SD) 

(49  -  82) 

(30-73) 

(49  -  60) 

1983  Canadian  Air 

98 

c 

2.4 

2.5 

2.6 

Sur-D  02  (DC8ASURW) 

O 

(1-3) 

(2-3) 

(2-3) 

With  input  from  NAVSEA  00C  and  the  Fleet,  we  evaluated  the  DCS  risk  associated  with 
possible  changes  to  the  present  Sur-D  O2  decompression  tables.  Because  no  single 
model  appeared  “best,”  we  used  both  the  ASYM  and  JAP98  models  to  evaluate  three 
possible  alternative  Sur-D  O2  procedures,  each  having  more  flexibility  than  current 
procedures  but  varying  in  its  amount  of  DCS  risk.  (We  have  been  told  that  until  the  final 
pull  to  the  surface  for  chamber  decompression,  the  flexibility  of  having  Sur-D  O2 
schedules  match  the  in-water  air  schedules  for  the  same  range  of  depths  and  bottom 
times  would  benefit  dive  supervisors.)  As  the  calibration  data  has  been  expanded  with 
this  modeling  effort,  we  expected  the  risk  predictions  from  the  ASYM  models  to  differ 
from  the  predictions  using  JAP98.  However,  we  gained  added  confidence  when  and 
where  there  was  agreement  between  the  models. 

As  a  first  step,  we  wanted  to  evaluate  how  risky  current  Sur-D  O2  procedures  are.  The 
first  dive  we  considered  was  to  120  fsw  for  60  min,  followed  by  decompression  using 
the  current  120/60  schedule  (Table  6).  The  Fleet  (i.e.,  U.S.  Navy  master  divers) 
indicated  that  this  dive  profile  was  uncomfortably  risky  and  is  not  routinely  used.  As 
illustrated  in  Table  6,  the  predicted  risk  ranges  from  5%  to  7%,  which  agrees  with  the 
Fleet’s  perception  of  high  risk.  The  second  dive,  to  1 12  fsw  for  52  min,  demonstrates 
the  range  of  risk  possible  due  to  the  step  sizes  (time  and  depth)  of  the  current  Sur-D  O2 
Table;  this  range  demonstrates  that  the  risk  can  be  changed  by  >  2  %  while  diving  a 
correctly  chosen  schedule.  The  third  dive  “jumps”  to  the  120  fsw  for  90  min  schedule,  a 
dive  that  was  reportedly  made  during  the  TWA  800  recovery  operation  and  that  cuts  the 
risk  by  about  one-third.  With  regard  to  these  models,  the  2.5%  to  3%  risk  of  this 
schedule  might  be  viewed  as  an  estimate  of  the  degree  of  risk  acceptable  to  the  Fleet. 
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Table  6.  Predicted  DCS  risks  of  surface  decompression  schedules. 


Dives  (fsw/min), 
Decompression  Schedule 

Risk  Predicted  by  Model  (%) 

(95%  Confidence  Interval) 

JAP98 

ASYM 

Human/Rat 

ASYM 

Human-only 

120/60  on  120  /  60  Sur-D  02  Schedule 

7.0 

(5.8 -8.1) 

4.8 

(4.0 -5.5) 

5.2 

(4.5 -5.9) 

1 1 2  /  52  on  1 20  /  60  Sur-D  02  Schedule 

4.7 

(3.8 -5.7) 

3.7 

(3.0 -4.3) 

4.0 

(3.4 -4.5) 

1 1 2  /  52  on  1 20  /  90  Sur-D  02  Schedule 

3.2 

(2.2 -4.2) 

2.6 

(1.9-3. 3) 

2.5 

(1.9-3.0) 

We  next  wanted  to  examine  how  these  models  behaved  with  respect  to  different  types 
of  dive  profiles  before  we  proceeded  to  possible  new  Sur-D  O2  options.  Unfortunately, 
the  use  and  evaluation  of  probabilistic  decompression  models  can  quickly  produce  a 
collection  of  numbers  far  too  large  to  manage.  Thus,  we  adopted  a  standard  set  of 
hypothetical  profiles  to  characterize  model  performance  (Table  7).  Previous  studies 
have  shown  some  of  the  standard  profiles  included  here  to  be  particularly  "safe”  or 
particularly  "dangerous.”  We  include  brief  descriptions  of  real  dives  similar  to  the 
hypothetical  profiles  to  assess  the  models. 

The  first  two  profiles  in  Table  7  are  for  submarine  escape  profiles.  Use  of  the  British 
SEIE  equipment  has  been  uneventful  in  many  exposures  to  400  fsw  and  shallower. 
However,  deeper  trials  have  produced  some  cases  of  DCS.20  21 

The  60/60  with  “rapid”  ascent  is  a  common  benchmark.  Some  people  think  that  slowing 
the  ascent  to  take  an  additional  minute  or  two  would  significantly  reduce  DCS  risk,  but 
the  value  of  this  approach  remains  unsubstantiated.  Many  thought  that  an  additional  20 
min  at  60  fsw  would  be  devastating,  but  it  did  not  appear  so  when  actually  performed.22 

The  60-min  exposures  to  49  and  80  fsw  have  the  same  nitrogen  profile  as  60  fsw  on  air, 
but  the  oxygen  partial  pressure  varies,  a  procedure  previously  shown  to  have  negligible 
effects.23 

The  three  saturation  profiles  span  the  range  from  no  observed  DCS  (20  fsw)  to  a  non¬ 
trivial  observed  incidence  of  DCS  (30  fsw). 

The  initial  60  fsw  for  three  hours  schedule  was  problematic  and  required  substantially 
more  decompression  time  than  required  by  the  air  tables  to  become  acceptable.24  The 
next  several  profiles  represent  alternate  ways  to  decompress  from  the  same  exposure. 
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The  final  pair  is  surface  decompression  possibilities,  one  of  which  was  found  more 
satisfactory  than  another.1 

From  Table  7  we  conclude  that  the  estimates  are  not  fully  consistent  with  experience 
(See  Appendix  F  for  confidence  intervals).  As  with  the  prior  model  USN93,  all  three 
models  overestimate  the  risk  from  a  20  fsw  saturation  dive.  Overestimation  probably 
occurs  for  the  shortest  dives  as  well  as  for  rapid  submarine  escapes.  Within  the  middle 
time  range,  it  appears  that  ~  3%  is  "safe"  by  JAP98,  while  >  6%  is  probably  "unsafe". 
The  two  ASYM  models  behave  similarly:  they  also  seem  to  have  ~  3%  as  "safe"  and 
>  6%  as  “unsafe.”  Other  benchmarks  can  be  used  based  upon  the  predictions  for 
schedules  that  are  familiar  (Appendix  B). 

Table  7.  Standard  profiles  with  predictions  from  the  three  models. 


Profile  Description  (All  air  dives,  unless  noted) 

Risk  Predicted  b' 

Model  (%) 

JAP98 

ASYM 

Human/R 

at 

ASYM 

Human-orily 

Submarine  escape  from  400  fsw  in  SEIE 

5.0 

8.2 

8.6 

Submarine  escape  from  600  fsw  in  SEIE  (2.7  min) 

10.1 

13.1 

14.3 

60  fsw  for  60  min  ascent  60  fsw/min 

2.2 

2.5 

2.3 

60  fsw  for  60  min  ascent  20  fsw/min 

2.1 

2.2 

2.1 

60  fsw  for  80  min  ascent  60  fsw/min 

3.3 

3.5 

3.3 

60  fsw  for  80  min  ascent  20  fsw/min 

3.2 

3.1 

3.1 

49  fsw  for  60  min  on  10%  Oxygen  (60/60  equiv) 

2.3 

2.4 

2.2 

80  fsw  for  60  min  on  35%  Oxygen  (60/60  equiv) 

2.2 

2.8 

2.5 

20  fsw  saturation  ascent  60  fsw/min 

6.8 

5.8 

5.1 

25  fsw  saturation  ascent  60  fsw/min 

10.9 

12.3 

9.2 

30  fsw  saturation  ascent  60  fsw/min 

15.2 

21.9 

19.0 

60  fsw  for  1 80  min  by  USN  in-water  schedule 

9.3 

6.9 

7.0 

60  fsw  for  180  min  by  USN93  schedule 

6.4 

5.7 

5.9 

60  fsw  for  180  min  with  14  min  O2  at  20  fsw 

8.0 

6.6 

6.7 

60  fsw  for  180  min  with  14  min02  during  ascent 

8.0 

6.5 

6.7 

60  fsw  for  180  min  by  USN  Sur-D  O2  schedule  — 

1  min  SI 

5.6 

5.5 

5.4 

60  fsw  for  180  min  by  USN  Sur-D  02  schedule  — 

3.5  min  SI 

5.7 

5.5 

5.5 

60  fsw  for  180  min  by  USN  Sur-D  02  schedule  — 

6  min  SI 

5.7 

5.5 

5.6 

60  fsw  for  180  min  —  1  hour  Oxygen  at  20ft 

2.5 

4.2 

4.2 

TWA  800—  11 6/85  on  the  1 20  for  90  Sur-D  02  schedule 

6.7 

5.3 

5.5 

TWA  800  —  1 1 6/59  on  the  1 20  for  90  Sur-D  02  schedule 

3.2 

3.2 

3.2 
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We  have  no  reason  to  believe  that  one  of  these  models  is  better  than  the  others.  We 
proceeded  with  all  3  models,  knowing  that  if  the  predictions  diverged,  we  were 
extrapolating  beyond  the  applicable  depth  time  range  of  at  least  one  of  the  models,  we 
would  then  need  additional  tools. 

We  then  began  to  use  the  three  models  to  predict  the  risk  of  possible  changes  to 
current  U.S.  Navy  Sur-D  02  procedures.  One  of  the  issues  we  wished  to  address  was 
the  apparent  need  by  U.S.  Navy  divers  to  jump  Sur-D  02  decompression  schedules 
(e.g.,  to  use  a  120/90  decompression  schedule  following  a  120/60  dive)1.  At  the  same 
time,  we  wanted  to  add  the  flexibility  of  having  Sur-D  02  schedules  match  the  air 
schedule  until  the  final  pull  to  the  surface  for  chamber  decompression,  flexibility  that  will 
allow  divers  easily  to  switch  from  the  Navy's  standard  air  decompression  to  Sur-D  02 
without  needing  to  decide  before  starting  decompression,  as  they  now  must  do.  Lastly, 
we  wanted  to  expand  the  procedures  to  cover  the  same  depth  and  time  ranges  as  those 
for  in-water  decompression.  Currently,  the  Sur-D  02  tables  go  only  to  a  maximum  depth 
of  170  fsw,  while  the  air  decompression  tables  continue  to  190  fsw. 

As  a  first  step,  we  assessed  the  risks  of  schedules  in  both  the  present  air 
decompression  tables  and  the  present  Sur-D  02  tables  (Appendix  B).  We  then 
examined  three  possible  Sur-D  initiation  rules  from  standard  air  (STD  AIR)  schedules  as 
demonstrations  for  selecting  new  surface  decompression  schedules.  This  approach 
may  allow  a  decision  about  adopting  a  given  procedure  for  Fleet  use  to  be  made  without 
requiring  a  costly  manned  dive  trial. 

The  three  Sur-D  initiation  rules  from  STD  AIR  schedules  we  chose  were  the  following: 

1)  4  0/m  in  5/ha  If.  Perform  all  the  in-water  decompression  required  by  the  air 
tables  until  the  40  fsw  stop  is  reached.  Then  perform  either  a  full  or  a  five- 
minute  stop  time  (whichever  is  greater),  and  pull  to  the  surface.  Return  the 
diver  to  pressure  in  a  chamber  within  five  minutes,  and  perform 
decompression  on  02  at  40  fsw  equivalent  for  half  of  the  remaining 
decompression  time  (30,  20,  and  10  fsw  stops).  Air  breaks  are  to  be  used, 
but  these  times  do  not  count  toward  the  02  decompression  time. 

2)  30/min5/half.  Perform  all  the  in-water  decompression  required  by  the  air 
tables  until  the  30  fsw  stop  is  reached.  Then  perform  either  a  full  or  a  five- 
minute  stop  time  (whichever  is  greater),  and  pull  to  the  surface.  Return  the 
diver  to  pressure  in  a  chamber  within  five  minutes,  and  perform 
decompression  on  02  at  40  fsw  equivalent  for  half  of  the  remaining 
decompression  time  (20  and  10  fsw  stops).  Air  breaks  are  to  be  used,  but 
these  times  do  not  count  toward  the  02  decompression  time. 

3)  40/min5/three-quarters.  Perform  all  the  in-water  decompression  required  by 
the  air  tables  until  the  40  fsw  stop  is  reached.  Then  perform  either  a  full  or  a 
five-minute  stop  time  (whichever  is  greater),  and  pull  to  the  surface.  Return 
the  diver  to  pressure  in  a  chamber  within  five  minutes,  and  perform 
decompression  on  02  at  40  fsw  equivalent  for  three-quarters  of  the  remaining 
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decompression  time  (30,  20,  and  10  fsw  stops).  Air  breaks  are  to  be  used, 
but  these  times  do  not  count  toward  the  O2  decompression  time. 

The  criteria  to  be  compared  among  these  options  and  the  present  Sur-D  02  procedures 
include  predicted  risk  of  DCS  and  total  decompression  time  (TDT).  Other  factors  that 
may  be  important  are  the  total  in-water  decompression  time  and  the  amount  of  time 
spent  on  oxygen  because  of  concerns  about  oxygen  toxicity.  The  amount  of  in-water 
time  will  depend  on  the  depth  of  the  last  in-water  stop.  The  risk  of  oxygen  toxicity  will 
be  assessed  after  a  procedure  has  been  chosen;  this  risk  can  be  minimized  by  using  air 
breaks  during  the  chamber  decompression. 


Table  8.  Risk  comparison  of  surface  decompression  schedules. 


Dives  (fsw/min), 
Decompression  Schedule 

Risk  Predicted  by  Model  (%) 

(95%  Confidence  Interval) 

JAP98 

ASYM 

Human/Rat 

ASYM 

Human-only 

120  /  60  on  120  /  60  Sur-D  02  Schedule 

7.0 

(5.8  -  8.1 ) 

4.8 

(4.0 -5.5) 

5.2 

(4.5-  5.9) 

1 12  /  52  on  120  /  60  Sur-D  02  Schedule 

4.7 

(3.8 -5.7) 

3.7 

(3.0 -4.3) 

4.0 

(34-4.5) 

1 1 2  /  52  on  1 20  /  90  Sur-D  02  Schedule 

3.2 

(2.2 -4.2) 

2.6 

(1.9 -3.3) 

2.5 

(1.9-  3.0) 

1 1 2  /  52  on  Option  1 :  40/min5/half 

4.6 

(3.6 -5.5) 

3.4 

(2.4  -4.4) 

3.4 

(2.8 -4.0) 

1 12  /  52  on  Option  2:  30/min5/half 

4.5 

(3.5 -5.5) 

3.2 

(2.1  -4.3) 

3.3 

(2.7-  3.9) 

1 1 2  /  52  on  Option  3:  40/min5/three-quarters 

3.5 

(2.5 -4.4) 

3.0 

(2.2 -3.9) 

3.0 

(24-3.6) 

All  three  of  the  Sur-D  02  options  we  examined  have  smaller  TDT s  than  those  of  the 
present  air  tables  (which  were  deliberate  choices;  see  Appendix  C).  However,  option  3 
(surface  from  40  fsw,  3/4  of  the  remaining  decompression  time)  is  only  slightly  shorter 
than  the  present  Sur-D  02  schedules,  with  its  advantages  being  shorter  in-water  time, 
and  decreased  risk  of  DCS  (Appendices  D  and  E).  With  a  maximum  predicted  risk  of 
6.5%,  this  option  has  an  appreciably  lower  maximum  risk  than  those  of  the  other  two 
options  and  the  present  air  or  Sur-D  02  tables:  9%  for  the  present  air,  and  1 1%  for  the 
Sur-D  02  tables.  Thus,  option  3  is  the  most  conservative  choice  for  risk  of  DCS. 

Options  1  and  2  have  similar  TDTs,  with  predicted  risks  slightly  lower  for  option  2 
(surface  from  30  fsw)  than  option  1 .  The  big  difference  between  these  two  choices  is  in 
the  amount  of  in-water  decompression  time.  With  both  of  these  options,  the  TDTs  are 
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up  to  10%  longer  and  the  risks  of  DCS  are  slightly  less  than  those  of  the  present  Sur-D 
02  procedures. 

We  can  now  examine  these  three  new  options  in  relation  to  the  1 12  fsw  for  52  min 
profile  (Table  8).  In  this  case,  we  can  see  that  all  three  options  have  predicted  risks 
lower  than  the  present  schedule,  but  not  as  low  as  jumping  to  a  90-minute  schedule. 

The  first  two  options  have  very  similar  risk  predictions,  with  their  differences  being  in  the 
amount  of  time  spent  in  the  water.  Using  the  ASYM  models,  we  can  also  see  that  the 
difference  between  the  options  is  minimal,  although  the  TDT  varies  from  42  to  60 
minutes  for  these  two  options.  This  18-minute  difference  in  decompression,  even  while 
on  oxygen,  affects  the  predicted  outcome  of  the  dives  only  by  <  0.5%  based  on  the 
ASYM  models;  the  difference  is  greater  for  JAP98. 

Using  the  TDT  information  plus  the  risk  predictions  compared  to  the  present  air  and 
Sur-D  02  tables  in  Appendices  D  and  E,  an  informed  risk  analysis  decision  can  be 
made  about  the  options  we  are  presenting. 

CONCLUSIONS 

There  was  no  improvement  to  the  predictions  of  the  ASYM  model  by  using  a  data  set 
combining  rat  and  human  data,  over  using  the  human  data  alone  when  calibrating  the 
model. 

The  value  of  these  models  will  be  in  their  future  use  as  additional  risk  prediction  tools 
that  help  in  determining  when  predictions  are  extrapolations. 
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Appendix  A 


Predicted  DCS  Risks  of  Data  Sets  in  Reference  8 
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_  Predicted  DCS  Risks  of  Data  Sets  in  Reference  8 _ 

_ _ Predictions _ _ 

Data  Set _ j  N  1  Observed  JAP98  1  ASYM  -  H/R  ]  ASYM  -  H 

Predictions  Confidence  Limits  Predictions  Confidence  Limits  Predictions  Confidence  Limits 
Single  Air _  _  _ _ 
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Appendix  B 


Predicted  DCS  Risks  of  the  Present  U.S.  Navy  Decompression 

Tables 

a)  Air  Tables 

b)  Sur-D  02  Tables 

Using  the  JAP98,  ASYM  Human/Rat,  (ASYM-H/R)  &  ASYM 
Human-only  (ASYM-H)  Models 
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Appendix  C 


Comparison  of  Decompression  Times  of  Proposed  Options 
a)  with  present  Air  Tables  (USN  Air) 
b)  with  present  Sur-D  02  Tables  (USN  Sur-D  02) 
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Appendix  D 


Risk  Comparison  of  Proposed  Options  with  Present  Air  Tables 
Using  the  JAP98,  ASYM  Human/Rat,  (ASYM-H/R)  &  ASYM 


Human-only  (ASYM-H)  Models 
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Risk  Comparison  of  Proposed  Sur-D  02  Options  Using  JAP98 


r\i 


(%)  *S!H  papipsJd 


Schedule 


Risk  Comparison  of  Proposed  Sur-D  02  Options  Using  ASYM-H/R 


CO 


(%)  >|S|d  papipaJd 


Schedule 


Risk  Comparison  of  Proposed  Sur-D  02  Options  Using  ASYM-H 


(%)  >iS!d  papjpaJd 


Schedule 


Appendix  E 


Risk  Comparison  of  Proposed  Options 
with  Present  Sur-D  O2  Tables 

Using  the  JAP98,  ASYM  Human/Rat,  (ASYM-H/R)  &  ASYM 
Human-only  (ASYM-H)  Models 
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Risk  Comparison  of  Proposed  Sur-D  02  Options  Using  JAP98 
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Risk  Comparison  of  Proposed  Sur-D  02  Options  Using  ASYM-H/R 
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Risk  Comparison  of  Proposed  Sur-D  02  Options  Using  ASYM-H 
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Appendix  F 


Standard  Profiles  with  95%  Confidence  Intervals 
for  the  Predictions  of  the  Three  Models 
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Standard  Profiles  with  Confidence  Intervals  for  the  Predictions  of  the  Three 

Models. 
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JAP98 

ASYM 

Human/Rat 

ASYM 

Human-only 

Submarine  escape  from  400  fsw  in  SEIE 

2.7 -7.1 

4.5  -  11.8 

4.3-12.7 

Submarine  escape  from  600  fsw  in  SEIE  (2.7  min) 

6.9-13.1 

7.6  -  18.3 

7.2-20.8 

60  fsw  for  60  min  ascent  60  fsw/min 

1.5 -2.9 

2.1  -2.9 

1.9 -2.7 

60  fsw  for  60  min  ascent  20  fsw/min 

1.4 -2.8 

1.7  -2.6 

1.6 -2.5 

60  fsw  for  80  min  ascent  60  fsw/min 

2. 5-4.1 

2.9  -4.0 

2.8 -3.7 

60  fsw  for  80  min  ascent  20  fsw/min 

2.4 -4.0 

2.6  -3.6 

2.6 -3.6 

49  fsw  for  60  min  on  10%  Oxygen  (60/60  equiv) 

1.6 -3.0 

1.9  -2.9 

1.8 -2.6 

80  fsw  for  60  min  on  35%  Oxygen  (60/60  equiv) 

1.0-  3.3 

2.2 -3.4 

2.1  -3.0 

20  fsw  saturation  ascent  60  fsw/min 

5.4 -8.2 

4.4  -7.3 

4.2 -5.9 

25  fsw  saturation  ascent  60  fsw/min 

8.7-13.0 

9.0-15.5  5.5-12.8 

30  fsw  saturation  ascent  60  fsw/min 

12.0-18.2 

16.3-27.0  12.4-25.1 

60  fsw  for  1 80  min  by  USN  in-water  schedule 

7.4-11.1 

5.9  -7.9 

5.9 -8.0 

60  fsw  for  180  min  by  USN93  schedule 

5.0 -7.7 

4.9  -6.6 

5.0 -6.8 

60  fsw  for  180  min  with  14  min  O2  at  20  fsw 

6.5 -9.6 

5.5 -7.7 

5.7 -7.6 

60  fsw  for  180  min  with  14  min02  during  ascent 

6.4 -9.6 

5.4  -7.5 

5.7 -7.6 

60  fsw  for  1 80  min  by  USN  Sur-D  O2  schedule  — 

1  min  SI 

4.2-7. 1 

4.6  -6.3 

4.6 -6.1 

60  fsw  for  180  min  by  USN  Sur-D  02  schedule  — 

3.5  min  SI 

4.3 -7.1 

4.6  -6.4 

47-6.4 

60  fsw  for  180  min  by  USN  Sur-D  O2  schedule  — 

6  rriin  SI 

4.3 -7.1 

4.6  -6.4 

47-6.5 

60  fsw  for  180  min  —  1  hour  Oxygen  at  20ft 

1.1  -3.9 

3.2  -5.3 

3.6-47 

TWA  800  —  11 6/85  on  the  1 20  for  90 

Sur-D  02  schedule 

5.2 -8.1 

4.3 -6.2 

4.6 -6.3 

TWA  800—  11 6/59  on  the  1 20  for  90 

Sur-D  O2  schedule 

2.1  -4.2 

2.5 -4.0 

CO 

1  00 

1 

CD 

CNJ 

1 
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